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' We analyze the possibility of universality violation in diagonal leptonic decays of 

the Z boson, in the context of interfamily "see-saw" models. In a minimal extension of 
the Standard Model with right-handed neutrino fields, we find that universality-breaking 
effects increase quadratically with the heavy Majorana neutrino mass and may be observed 
in the running LEP experiments. 
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Heavy Majorana neutrinos with masses of few TeV have entered the domain of cos- 
mology and astrophysics, as possible candidates to account for the net baryon number 
of the universe through lepton (L)-number violating processes On the other hand, 
if such heavy neutral leptons are realized in nature, their existence may be discovered 
through their production and the L-violating decay at present or future collider ma- 
chines [0] . Another place of looking for new physics originating from heavy neutrinos is 
the rare leptonic decays of and Z particle Since such decays are forbidden 
in the minimal Standard Model {SM), they constitute an interesting framework to con- 
strain new physics beyond the SM. Such rare decays |^ have recently been analyzed 
in a "see-saw" -type model with intergenerational mixings 0. An interesting aspect 
of this minimal scenario is that the Appelquist-Carazzone theorem is not operative 
and vertex-correction diagrams with intermediate heavy Majorana neutrinos (Ni) show 
a quadratic mass dependence, i.e. m^./M^. This mass dependence is a common feature 
for all theories based on the spontaneous-symmetry breaking mechanism. For example, in 
the Feynman-'t Hooft gauge this mass dependence can be seen to arise from the coupling 
of the would-be charged Goldstone bosons to heavy fermions. Similar non-decoupling ef- 
fects in the quark sector have been extensively studied in the past for the flavor changing 
decays Z — > 6s and the diagonal Z bb [^. 

In this note we study universality-breaking effects induced by heavy Majorana neu- 
trinos in leptonic Z decays. Actually, we will analytically calculate the following physical 
observable: 

^ r(Z^r+r-) - r(Z ^ /+/-,/ = e or /x) 

T{Z T+T-) + T{Z 1+1-, / = e or ^) ■ ^ ' 

Ubr is a measure of universality violation in the leptonic sector provided lepton masses can 
be neglected and is experimentally constrained to be [rT| \Ubr\ < 1-5 10^^. The model we 



are concerned with extends the SM by one right-handed neutrino field for each family. 
The renormalizable form of all relevant interactions and details of the notation we will 



use here can be found in ^. 

In the SM the universahty- violating parameter Ubr has a value different from zero 
due to r lepton mass effects. However, pure phase-space {PS) corrections turn out to 
be rather small, i.e. \ul!^^\ ^ 3 [1 + (1 - 4 sin^ ^vk)^]~^ m^/Mf ~ 1.1 10^1 This 
standard source for a non-zero value of U^r is beyond the sensitivity of the present LEP 
experiments. Since one can expect to analyze about 10^ leptonic decays of the Z boson 
at LEP per year, one may reach an accuracy for \U},r \ at the level of 3. 10~^. Since we are 

\ IPS')] 

interested in values of \Uhr\ much larger than \Ul^ |, we will neglect lepton-mass effects 
in the calculation of quantum corrections. 

Decomposing now the transition matrix element T{Z —>■ II) into two parts (with the 
superscripts (0), (1) denoting the electroweak loop order), i.e. T(Z II) = T^^\Z 
ll) + T^^\Z 11) , with 

rW(Z^/r) = ^^4n,7,(l-4sin2^H'-75K, 
4 cos 6w 

and defining AT^^) = T^^\Z r+r") - T^^\Z l+r, / ^ r) , we find that 

Ar(i) = '^^"^ 4 ^^7.(1 - l,)v! AB,,Fz{\, A,) . (2) 
8tt cos Ow 

Then, the universality-breaking parameter Ubr takes the simple form 

- fPfi)[^ - V l + (l-4sin2M^ ' ^'^ ' ^ ' 



where 
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m: 



= T7T' ' ^^ij = B^iB*,: - BuBi, , I ^ t . (4) 

In Eq. (4), Bu is a Cabbibo-Kobayashi-Maskawa-type hq x 2nG matrix appearing in the 
leptonic charged-current interaction, and indicates the masses of all neutral leptons rii 
in our minimal scenario (i.e. i = 1, 2, . . . , 2nGi with no denoting the number of families). 
The function Fz{Xi, Xj) originates from the one-loop graphs depicted in Fig. 1 and contains 



all the non-decoupling physics mediated by heavy Majorana neutrinos. Its explicit form 
will be discussed below. 

Before proceeding to give the analytical form of Fz-, we remark that the universality- 
violating parameter Ubr does not involve any infra-red {IB) singularities, as they do not 
depend on neutrino masses. Therefore, soft-photon emmision graphs need not to be 
considered here. Also, the Z — I — I vertex should be renormalized at the one loop level. 
In fact, a large number of renormalization constants should not appear in the universality- 



breaking parameter ^7^.. If we adopt the on-shell renormalization scheme |T2| , |T^, where 
the input renormalization parameters are the electric charge e. My/-, Mz, the Higgs mass 



Mh and all fermion masses contained in the model, the counterterm Lagrangian C 



c 

int 



relevant for the renormalization of the Z — I — I vertex is then given by |n 



c 

int 



le Z'il hi^ 2^ + 



le 



4:SwCw 



5e 1 
1 + — + - 

e 



■~ 2c^ 112 
— bp + bZzz 



+ 5Z[ 



Cw 



5e 1 
1 + — + 

e 



9^2 



5p 



bzT 



75 



(5) 



2 ^2 

where bp = SM^/Mz — bMy^/M^^, cw = cos^^^ = Mw/Mz and sw = sin6'H^. Due to 



G/M-type cancellation the only non- vanishing contribution to the function Fz comes 
from the wave-function renormalization constants of the left- and right-handed leptons, 
i.e. bZ^i^ and bZ^^. In fact, one has to calculate the difference of the self-energy derivatives 
given by A = bZ}^ — bZ^j^^ = A dT,(^) / d^\j^_^ m ^o' The corresponding constant for 
the right-handed leptons AZjj, vanishes in the limit mi, — > 0. It is easy to see that only 
the neutrino-mass dependent self-energy graphs mediated by and are of interest 
here. The individual contributions to Fz{Xi,Xj) arising from the diagrams l(a)-l(f) and 
those from the wave-function renormalization constant AZ^ are given by 



Cij ( L2{\i, \j 



A; 



(6) 
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CijXiXjKi{Xi, Xj) + C*jyXiXj[ 2^uv 



- + XzK{\,Xj: 



6ij 



XzhXi) + 34Li(A,) 



(7) 



(9) 



if) 



where 



S,, ^ A,/(A,) 

Cw 



-Ml-2sl,)X, C 



uv 



3 
2 



1 - A,- 



(A^ + 2)A^ In A, 
(1 - A.)^ 



no 



hi ' 



k=l 



Ml 



C 



uv 



7£; + ln47r — In 



(10) 
(11) 



(12) 



The functions /, /, Li, Ki, K2, K and L2 involved in Eqs. (6)-(ll) are given in Ap- 
pendix A. It is straightforward to see that the ultraviolet {UV) divergences (i.e. Cuv) 
cancel in the summation of all Fz terms. To be precise, the UV pole in F^'' cancels 
against the UV one of the wave-function renormalization ^ and the UV constant in 



2nG 

Eq. (7) vanishes due to the identity y]: X) ^n.BiiC*^ 

1=1 



. 



For definiteness, we will consider a interfamily-mixing model with two families 
only. Employing relations between the mixing matrix Bn and heavy Majorana neutrino 
masses p together with Eqs. (3) and (6)-(ll), we arrive at the simple result 



\u 



br 



<^w (sl )^ 



in 



-1/2)2 



X 



1 + llna; - (1 + 2x'/^) 

2 1 



X 



(13) 



where x = 'm%^/m'j^^. We also assume N2 to be heavier than Ni, i.e. a; > 1. In Eq. (13) 
s'2 is the usual neutrino-mixing angle between heavy Majorana neutrinos and the charged 



lepton / and is generally defined as (-Sl)^ — S IBin]"^ ■ The neutrino-mixing angle 
(s^^)^ turns out to be severely constrained by the recent LEP data on r decays. In 
fact, a global analysis allowing mixing of exotic particles gives an upper bound of about 
7. 10^^ ||T^. The e- and /x-family is much more constrained, i.e. {s^i^^^Y < 0.01. Due 
to this fact we have to deal with universality-breaking effects in the heaviest lepton fam- 
ily only. It is important to notice that the non- decoupling terms (i.e. proportional to 
m^/M^) come from the "seemingly" suppressed (s^^)^ terms. Table 1 shows the dra- 
matic non- decoupling behaviour of the loop function Fz- For comparison, we also show 
the corresponding values for a calculation where terms proportional to {s'^Y have been 
neglected. In our numerical estimates we have assumed that there is no large mass dif- 
ference between the two heavy Majorana neutrinos, i.e. x ~ 1. We have varied the heavy 
neutrino mass uin- (~ mAr^ ~ mATj) up to its perturbative unitarity bound. Such an upper 
bound may be imposed by requiring that the total width of Ni and A^2, denoted by T^i, 
satisfies the condition T^JniN. < 1/2 f^. This leads to the constraint on the mass of the 
lightest heavy Majorana neutrino A^^i 

""^^ - a^is^y ■ ^^^^ 

Taking the above upper bound into account, we find that the universality- violating pa- 
rameter \Ubr\ can be up to 10 times larger than the naive value obtained by considering 
only terms proportional to the mixing (s^^)^. This enhancement factor (i.e. the crucial 
m%/M^r dependence) results from the coupling of the would-be charged Goldstone boson 
X"^ to the heavy Majorana neutrinos in the diagram 1(b). If we assume very large mass 
differences between the heavy neutrinos, \Ubr\ smoothly decreases to negligibly small val- 
ues. The reason is that one effectively recovers the one generation "see-saw" model in 
such a case, i.e. m^Vi ^ as x — > oo in Eq. (14). In Fig. (2) we present exclusion plots 
for LEP experiments. We see, for example, that possible universality-breaking effects of 
the order of 10~^ can easily be understood within our minimal model. 
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Non-5'M contributions, due to heavy neutrinos, could also be observed by compar- 
ing the prediction of the SM for the leptonic rates T{Z — > l^l^) with those epxected when 
one includes heavy neutral lepton effects. Besides the universality-violating phenomena 
discussed in this paper, one could have non-decoupling physics introduced by these heavy 
neutral leptons which can be constrained by analyzing the oblique electroweak parameters 



5, T, f/ or ei, £2, ^3 as defined in [|T^. These contributions occur both through vacuum 
polarization terms and in vertex corrections universal for all charged lepton flavors. How- 
ever, a first estimate for the 5p parameter |]19[ (i.e. 5p = aemT) provides a weaker bound 
than that obtained by Eq. (14). In particular, the dominant non-S'M contribution (de- 
noted below as 6p*) to 6p, apart from heavy-top and Higgs-particle effects, comes from 
the I'iNj, NiNj intermediate states of the Z self-energy graph. In this way, one gets 

(15) 

We can readily compare Eq. (15) with constraints on the masses of the heavy neutrinos 
that are derived on the basis of perturbative unitarity. For example, if (s^^)^ = 0.1, the 
maximal value that 6p* can take is 0.8 10~^, which is still in accordance with phenomeno- 



logical contraints of a possible mass shift of the W boson However, to complete the 
analysis a global consideration of all electroweak oblique corrections is required. 

In conclusion, we have explicitly shown that heavy neutral leptons introduce a 
quadratic mass dependence (i.e. awfn%/M^r) in the leptonic vertex function Z — I — I. 
This situation is not peculiar for the minimal model considered in this work, but a gen- 
eral feature for all theories based on the spontaneous symmetry-breaking mechanism, e.g. 
similar effects will be present in the model described in [Q. In general, we have found that 
the mass of possible non-decoupling neutral particles can be constrained by the already 
existing or future LEP data, as a function of their mixing to the ordinary charged leptons. 
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A The one-loop integrals 

It is useful first to define the functions 5i(Ai) and B2{Xi, Xj) which are given by the 
following expressions: 

B^(Xi) = (1 - y)Xi + y[l - Xz yx(l - x)] , (Al) 

B2{Xi,Xj) = l-y + y[xXi + {l-x)Xj- Xzyx{l-x)] , (A2) 

where x and y are Feynman parameters. The loop functions /, /, Li, Ki, K2, K and L2 
can then be written in terms of the following integrals: 



Li(A,) = J dxdyy In B,{Xi) , (A5) 

- J 5,(A„A,) ' ^^^^ 
L2(Ai, Xj) = J dxdy y In B2{Xi, Xj) , (A9) 

where the integration variables x and y are constrained to the interval [0, 1] . The above 
integrals can be best performed numerically. 
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Figure and Table Captions 



Fig. 1: One-loop irreducible vertex graphs contributing to the non- universality pa- 
rameter Ubr in the Feynman-'t Hooft gauge. 

Fig. 2: Exclusion plots for LEP experiments. The areas lying to the right of the 
curves are excluded due to the following conditions: (i) The validity of 
perturbative unitarity (sohd hne) - see also text, (ii) |C4r| < 1- 10~^ (dashed 
line), (iii) \Ui,r\ < 7. lO'^ (dash-dotted line), (iv) \Ubr\ < 3. lO'^ (dotted 
line) . 

Tab. 1: Numerical estimates for the universality- violating quantity \Ubr\ within the 
perturbatively allowed parameter space. The values in parentheses are 
obtained by neglecting the "seemingly" suppressed terms proportional to 



11 



Table 1 
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